Dipolar depletion effect on the differential capacitance of carbon based materials 
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The remarkably low experimental values of the capacitance data of carbon based materials in 
contact with water solvent needs to be explained from a microscopic theory in order to optimize 
the efficiency of these materials. We show that this experimental result can be explained by the 
dielectric screening deficiency of the electrostatic potential, which in turn results from the interfacial 
solvent depletion effect driven by image dipole interactions. We show this by deriving from the 
microscopic system Hamiltonian a non-mean-field dipolar Poisson-Boltzmann equation. This can 
account for the interaction of solvent molecules with their electrostatic image resulting from the 
dielectric discontinuity between the solvent medium and the substrate. The predictions of the 
extended dipolar Poisson-Boltzmann equation for the differential capacitance are compared with 
experimental data and good agreement is found without any fitting parameters. 

PACS numbers: 03.50.De,05.70.Np,87.16.D- 
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I. INTRODUCTION 

New generation supercapacitors are used for a broad 
range of applications in nanoscopic scale technologies. In 
water purification technology, capacitive desalination is 
an efficient candidate that might replace the current lead- 
ing technics such as reverse osmosis, a membrane based 
purification method known to suffer from the membrane 
fouling phenomenon [Tj. Supercapacitors are also used 
as low cost and long life energy storage devices with con- 
siderably higher energy densities than conventional elec- 
trolytic capacitors [2|. A through understanding of the 
double layer structure of these devices is thus necessary 
to optimize their efficiency. 

The understanding of the double layer structure was 
limited for several decades to the Gouy-Chapman-Stern 
model [S]- This model was later completed by consider- 
ing additional effects specific to electrolyte solutions, to 
name but a few, the steric layer associated with the size of 
solvent molecules as well as the dipolar alignment close to 
the interface [J] , non-local effects in electrolytes at metal- 
lic interfaces [5] , ionic crowding [6] and overscreening [7] . 

Supercapacitors are commonly fabricated from carbon 
based materials with a dielectric permittivity e m sw 2 — 5 
much lower than the permittivity of the water solvent 
s w = 78. The polarization of the interface resulting from 
this large dielectric discontinuity can drastically change 
the physics of the double layer. Image dipole interac- 
tions were considered in Ref. [5] for a metallic interface. 
However, the work accounted exclusively for the effect 
of image interactions on the dipolar orientation without 
considering their role on the interfacial dipole density. 
Furthermore, it was recently shown in Ref. [5] that the 



Gouy-Chapman (GC) capacitance largely overestimates 
the experimental data obtained for carbon based mate- 
rials. The failure of the GC capacitance was explained 
by the unability of the Poisson Boltzmann formalism to 
account for non-local dielectric effects. 

In order to gain insight about the contribution of sur- 
face polarization effects on the capacitance of low dielec- 
tric substrates, we introduce in this work a first micro- 
scopic modeling of solvent molecules beyond the MF level 
approximation. Namely, we derive an extended dipolar 
PB (EDPB) equation that can self-consistently take into 
account the interfacial solvent depletion. This depletion 
results from the interaction of solvent molecules (modeled 
as dipoles) with their electrostatic image, an effect absent 
in the mean-field level DPB equation [TTJ[T2] . The predic- 
tion of the EDPB equation is shown to agree well with ex- 
perimental data for the differential capacitance of carbon 
based materials. However, it is also shown that the DPB 
formalism yields the same result as the PB equation, that 
is, it overestimates the experimental data by one order 
of magnitude. These observations strongly suggests that 
the dielectric discontinuity between the substrate and the 
solvent can solely explain the observed low values of the 
differential capacitance of carbon based materials, unlike 
the conclusion of Ref. [9. where it was argued that the 
surface polarity does not play a major role in the differen- 
tial capacitance. Our results are also in agreement with 
the experimental work in Ref. |10j , where the surface hy- 
drophobicity was actually shown to strongly reduce the 
capacitance of carbone nanotubes. 



II. EXTENDED DIPOLAR POISSON 
BOLTZMANN (EDPB) EQUATION 
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We will present in this part the derivation of an ex- 
tended dipolar Poisson-Boltzmann formalism. The field- 
theoretic partition function of ions immersed in a dipo- 



lar liquid was derived in Ref. QT] as a functional integral 
over a fluctuating electrostatic potential <p(r) in the form 
Z = J 'T>4> e~ H w\ where the Hamiltonian functional is 
given by 
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The first integral term of the Hamiltonian (TTJ) is com- 
posed of the Maxwell tensor associated with a freely prop- 
agating electric field V^(r) in the air, and a second part 
that couples the corresponding potential 0(r) to a fixed 
surface charge distribution c(r). The second and third 
integrals respectively account for the presence of solvent 
molecules (point dipoles) and ions of different species de- 
noted by the index i. Moreover, r = (x, y, z) is the config- 
urational space and fl = (9, ip) stands for the solid angle 
characterizing the orientation of solvent molecules, with 
9 the angle between the dipole and the z axis. We note 
that the external wall potential V w (r) in Eq. (fTl) restricts 
the space accessible to the particles, and in the case of the 
single dielectric interface located at z = 0, it is of the form 
V w (z < 0) = oo and V w (z > 0) = 0. Furthermore, Xd and 
Xi are respectively the fugacity of dipoles and ions, p the 
dipole moment vector, and Qi stands for the valency of 
ions for the species i. The heterogeneous Bjerrum length 
is defined as ^s(r) = e 2 / [47re(r)/csT], where e is the ele- 
mentary charge, T = 300 K is the ambient temperature, 
and e(r) = £ 9(z)+e m 9(— z) is the dielectric permittivity 
of the medium in the absence of solvent molecules for the 
same single planar interface geometry. More precisely, Eq 
and e m denote respectively the dielectric permittivity of 
the air (the subspace in z > 0) and the low dielectric 
subbstrate located at z < 0. From now on, the dielec- 
tric permittivities will be expressed in units of £o- We 
also note that the Bjerrum length in the air is £b ~ 55 
nm. Finally, the self energy of ions and polar molecules 
that are substracted from the potential and electrostatic 
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field respectively read Ei = ^v c {r — r')| r=r / and Ed — 
|(P" V r )(p- V r ')«c( r — r/ )i where the Coulomb operator 
in the air is defined as u* _1 (r,r / ) = - hj ^ AS(r - r'). 

In this work, we aim at investigating the model of 
Eq. (fTl) beyond the MF approximation where surface po- 
larization effects are absent [TTJ [T^] . One way to progress 
consists in opting for a variational minimization proce- 
dure that aims at finding the upper boundary for the 
dimensionless Grand potential of the system Q = — In Z 
by minimizing the variational Grand potential defined as 
Cl v = fio + (H — Hq) , where the reference Hamiltonian 
is a Gaussian functional of the form 
I 



Hn = 



[<P(v) - ^ (r)] Vo\r, r') [>(r') - i^(r')] 



(2) 



Furthermore, <f>o(z) is a variational external potential and 
the electrostatic trial kernel is chosen in the same form 
as in Refs. [T3HTS], 

VoH*y) = ^£ [-V(e v (r)V)+E v (ryl(r)} 6(r-r>), 

. (3) 
where the piecewise variational dielectric permittivity is 

defined as £„(r) = e w 6(z)+e m 9(— z) and the trial screen- 
ing length is given by K c (r) = k c 9(z). After performing 
the functional integrals over 4>{v), one gets 
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where the gaussian contribution reads Oo = 
-ln/P0 e - ff oM. We also defined above the local 
ion density 

(U(t) = A l e^- y "'( r )e- ? ^ o(r) -^ , ' o(r ' r) (5) 

and the local density of dipoles with orientation fi 
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By taking the derivative of the variational Grand po- 
tential Eq. (HI) with respect to n c and e v , one gets 
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^^BPoPbd- These two relations respectively introduce 
the Debye-Huckel screening parameter and the Debye- 
Langevin form for the bulk dielectric permittivity of the 
water medium e w . The additional variational equation 
for 4>o(z), i.e. 8Q v /8(po(r) = yields 

9 e{z)^^-+A-Kl B u{z)+A-Ki B Y d Pi{z)q i = 0, (7) 



dz 



dz 



where we took into account the translational symmetry 
of the electrostatic potential within the (x,y) plane. We 
note that the variational minimization left us in Eq. pj 
with a spatially varying dielectric permittivity of the 
form 



e(z) = 1 
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p d (z,ft)p z , 
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where p z = po cos 9 stands for the component of the 
dipolar moment vector p in the z direction. We will 
call Eq. (IT]) the Extended Dipolar Poisson Boltzmann 
(EDPB) equation. 

The fugacity of dipoles and ions can be related to their 
bulk density in the limit z — > oo of the equations ([5]) 



and ([6]). By injecting the obtained relations for the fu- 
gacities with the inverse of the kernel Eq. Efy [13] into 
Eqs. (O)]) and ([6]), the local density functions take the 
form 



Mz)=P b ^- V - {z) e- q ' Mz) ~ VJz) (9) 

p d (z, n) = pue-v-Me-vW)-™'™, (10) 

where we defined the following ionic and dipolar poten- 
tials, 
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and 

p c = \J n 2 . + k 2 . Carrying out the integral over 9 in 
Eq. (|8t with the dipole density Eq. (flOl) and the potential 
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(pc - V k )/(Pc + r)k), r] = e m /e v 



Eq. ffl2h, the local dielectric permittivity takes the form 



e(z) = 1 + ^-eBP 2 p db e~ v ^e- u ^J(z), 



(15) 



where we defined the function 

J(z) 



x e^T^r {Erf [*+(z)] 



8T 3 d /2 (z) 
3e~ T "^ smh\p <j>' (z)} 



2T d (z) 
and the potentials 
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The EDPB Eq. tffl has to be solved numerically with 
the ionic density profiles of Eq. Ety and the dielectric 
permittivity profile of Eq. (15 1. 

The second order differential equation tffl should be 
solved with the boundary conditions 4>q(z — > oo) = 
and 4>' (z —> + ) = 2e w /p,, where the second boundary 
condition valid over the parameter domain < e m < s w 
follows by integrating Eq. ^j\ in the close neighborhood 
of the interface, and noting that according to the dipo- 
lar potentials of Eqs. (TT3J) and (14 1, one has p d (0) = 
and e(0) = 1 on the boundary. We also note that in the 
limit where the potentials V c (z), U d {z), and T d {z) van- 
ish, EDPB equation. ^j\ reduces to the mean field DPB 
equation of Refs. [TTII12] . 

We finally note that the orientation averaged density 
of solvent molecules is obtained according to p d (r) = 



/ 4?Pd( r ) ^)- Evaluating the integral over 9, the solvent 
density takes the form 

Pd (z) = Pbd -^-e- v ^e-^e P M^ (18) 
x {Erf [*+(z)] + Erf [*_ (z)]} . 



III. NUMERICAL RESULTS 

We will investigate in this part the EDPB Eq. (17f for 
a symmetric electrolyte composed of two ion species of 
bulk densities Pb,i — pu and valency qi — q. All numer- 
ical results will be derived for monovalent ions (q = 1) 
in contact with a negatively charged planar surface, i.e. 
a(z) = —a s 5(z) with a s > 0. We also note that within 
the convention adopted in this article, the surface charge 
a a is expressed in units of the elementary charge e. More- 
over, the model parameters p d b and e w are taken the 
same as in Ref. [Sj. Namely, the bulk density of sol- 
vent molecules is p d \, — 50.8 M, which yields with the 
dipole moment po — 1 A the bulk dielectric permittivity 
e w = 71. 

The potential profile obtained from the numerical solu- 
tion of the EDPB Eq. ([7} for the parameters p bi = 0.1 M, 
e m = 1 and <j s = 0.01 nm~ 2 is reported in Fig. [TJa. One 
notices that the potential profile is composed of three re- 
gions, namely two successive layers close to the interface 
where <f>o(z) behaves as a linear function of z, and a third 
layer over which (f>o(z) exponentially decays. 

In order to understand the form of the potential profile, 
we illustrate in Fig. [lithe form of the dielectric permittiv- 



ity e(z) and the screening parameter n (z) 



n n e 



-V c {z) 

in the vanishing surface charge limit of Eq. ¥i\. It is 
seen that with increasing distance from the surface, the 
dielectric permittivity increases from the air permittiv- 
ity e(z) — 1 to the bulk permittivity e(z) = e w over a 
distance h ~ 2 A. We note that this dipolar exclusion ef- 
fect is mainly due to the interaction of solvent molecules 
with their electrostatic images. Then, one sees that this 
solvent depletion regime is followed by an ionic depletion 
regime of thickness d ~ 6 A, an effect known to originate 
from image charge interactions |13j . 

Inspired by the behaviour of e(z) and k(z) that re- 
sults from the interfacial depletion of solvent molecules 
and ions, we will introduce a restricted variational ansatz 
based on a piecewise trial solution for the electrostatic 
potential. We assume that cj)o(z) is the solution of 
Eq. (T7J) in the linear limit of weak surface charge, with 
e(z) = 6{h — z) +e w 9(z — h) and k(z) = K c 9(z — d), where 
the dipolar and ionic depletion lengths h and d are trial 
parameters that will be obtained from a numerical op- 
timization procedure of the Grand potential of Eq. Q. 
The solution of Eq. |7|) with the above piecewise dielec- 
tric permittivity and ion density profiles, and satisfying 
the continuity of the potential 4>${z) and the displace- 
ment field D(z) = £(z)4>' {z) at z — 0, z — h, and z — d, 
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FIG. 1: (Color online) (a) Electrostatic potential profile 
(<7 S = 0.01 nm -2 ) and (b) renormalized density and dielectric 
permittivity profiles for e w = 71 and pu — 0.1 M. The red 



line in (a) is from the restricted variational ansatz Eq. (191 



and the dashed black line corresponds to the solution of the 
EDPB equation. 



reads 

<Mz) = 
<Mz) = 



— [l + K c {d-h)} + —(z-h), 0<Z<h 
flK c [I 
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Numerical optimization yields h = 0.6 A and d = 2.3 
A. We note that due to the piecewise nature of the trial 
potential in Eqs. (19 1, these values correspond approxi- 
mately to half saturation densities. Figure [T] shows that 
the potential profile obtained from the numerical opti- 
mization agrees very well with the general form obtained 
from the numerical solution of the EDPB equation. In 



Eqs. ( 19 1, the first linear regime at < z < h corresponds 



to the solvent depletion layer resulting from image dipole 
interactions. This layer associated with dielectric screen- 
ing deficiency is responsible for an amplification of the 
PB prediction of the surface potential by a factor of five. 
The second and third intervals correspond respectively 
to the usual ionic depletion and diffuse layers [13 . The 
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FIG. 2: (Color online) (a) Differential capacitance against the 
bulk ion concentration for a s — 0, e m = 1, and e w — 71. 
The black circles are the experimental data, the red solid 
curve is the result of the EDPB equation, the red squares 
are from Eq. ( |21[ ), the dashed blue line is the MPB result, 
the dotted black line is the GC capacitance, and the black 
squares correspond to the prediction of the DPB equation. 
(b) The same plot as in (a) for various e m . The inset displays 
the evolution of the dipolar depletion length h (solid curve) 
and Id (dashed curve) as a function of e m for pti — 0.1 M. 



contribution of these layers to the differential capacitance 
will be investigated below. 

The differential capacitance of the double layer is de- 
fined as 



C,= 



k R T 



da* 



(20) 



where 4> s = 4>o{z = 0) is the surface potential. Fig. [2la 
compares the differential capacitance computed with 
Eq. (17| in the vanishing surface charge limit with ex- 
perimental data obtained for several types of monova- 
lent electrolytes at various concentrations (for details see 
Ref. [5] where the data were taken from). We also re- 
port in this figure the prediction of various formalisms. 
As stressed in Ref. [9], the PB result largely overesti- 
mates the experimental data. Furthermore, the result of 
the modified PB (MPB) equation (see Ref. [T3]) that can 
exclusively take into account the ionic depletion effect 



brings a very small correction to the PB result. How- 
ever, the EDPB result that additionally contains the sur- 
face depletion effect of solvent molecules exhibits a good 
agreement with the experimental data. We finally note 
that in the vanishing surface charge limit considered in 
this part, the DPB equation yields the same result as the 
PB one (see black squares in Fig.[2ja). 

The physics of the EDPB prediction for the capac- 
itance can be understood within the restricted self- 
consistent scheme of Eq. (19 1, where the differential ca- 
pacitance in Eq. ( 20 ) takes in the limit a s — > the simple 
form 



C,= 



Sqijhhr 



1 + K c (d — h) + £wKch 



(21) 



We note that the prediction of this equation reported 
in Fig. p] (red squares) fits very well the numerical re- 
sult of the EDPB equation. The inverse capacitance of 
Eq. ( pi) is composed of three parts. The first contri- 
bution from the diffuse layer is the inverse GC capaci- 
tance Cj} = (etoKc) -1 corresponding to the PB result 
in Fig. pi a. The second part C^ 1 = (d — h)/e w asso- 
ciated with the ionic depletion layer is shown to drop 
the differential capacitance to the MPB curve. Finally, 
the third contribution from the solvent depletion layer 
Cjz — h characterized by the dielectric screening defi- 
ciency brings the most important correction to the total 
capacitance by dropping the latter to the correct order 
of magnitude. 

In order to estimate the order of magnitude of the dipo- 
lar depletion length h, we will compute the asymptotic 
limit of e(z) far from the interface, where the dipolar po- 
tentials in Eqs. ( 13 1 and ( 14 1 become weak, by expanding 
Eq. (151 in Ud(zj~^Td(z), and po<fi' (z). Furthermore, we 
note that in the limit e m — 0, the dipolar potentials are 
given by the closed form expressions 



U d (z) — £ w pl 
T d {z) 



1 + 2hi r z 



-2k c z 



2 1 + 2k c z(1 + 2k c z) _ 2KaZ 
wPo 16z3 C 



(22) 
(23) 



Renormalizing all lengths by the length scale I4 — 
(IwPq/10) 1 ' 3 according to R c = nJd and z = z/ld, and 
taking into account that s w 3> 1, the asymptotic form of 
Eq. (Tl5|) far from the dielectric interface reads 



£ in 



1 _ l + 2^ + 3^z 2 /2 c „ 2BcZ - 



(24) 



We now note that for a bulk permittivity e w = 71 and 
ionic concentration p^ = 10 _1 M, one has n c ld = 8.10~ 2 . 
Hence, in the regime z ~ Id = 0.9 A, the terms in 
Eq. (24) that depend on the screening length become 



negligible. This simple calculation fixes Id as the charac- 
teristic length over which the local permittivity tends to 
its bulk value according to an inverse cubic power law, 
i.e. e(z)/e w wl- ^/z 3 - We note that an inverse cubic 
law for the dielectric permittivity profile was derived in 



Ref. [TB] in the strict limit of a single dipole (i.e. e w = 1) 
and without salt. In the dilute salt limit n c —¥ 0, one 
can actually extend the estimation of h to finite values of 
e m by noting that the dipolar potentials possed the close 
form expression Ud(z) = Td(z) — £ w PqAq/(16z 3 ), where 
Ao = (e w — e m )/(e w + e m ). Following the same steps as 
above, one obtains for the characteristic dipolar depletion 
length the more general expression Id = (Ao^p^/lO) 1 / 3 . 
Hence, for biological solvent concentrations and dilute 
electrolytes with bulk density pu < 0.1 M, the dielectric 
screening is mainly responsible for the decay of the image 
dipole interactions and solely determines the region over 
which a reduced dielectric permittivity is observed. For 
larger ion concentrations, Eq. (24) shows that the screen- 



ing of image dipole interactions by surrounding ions pos- 
itively adds to the dielectric screening of these forces. 

We display in the inset of Fig. [2]b the evolution of h 
as a function of e m together with Id, while the main plot 
shows Cd for various values of the membrane permittiv- 
ity from e m = 1 to e m = e w . One notices that within the 
range 1 < e m < 60, h exhibits a slow linear decrease with 
increasing e m while C'd remains within the same order of 
magnitude as the experimental capacitance data. How- 
ever, with an increase of e m from 60 to e w , the dipolar 
depiction length quickly drops to zero and consequently, 
Cd approaches the PB result. One also sees in the inset 
that although Id is slightly higher than the dipolar deple- 
tion length h, it can reproduce the correct trend of the 
latter as a function of e m . These observations suggest 
that image dipole interactions are mainly responsible for 
the low values of the experimental capacitance data in 
Fig. [2]a. We emphasize that this result is in agreement 
with the experimental observation of a strong reduction 
of the double layer capacitance with increasing surface 
hydrophobicity [10] . 

In addition to the dipolar depletion effect, the orienta- 
tion is also expected to play some role in the form of the 
interfacial dielectric permittivity profile. The measure 
of the dipolar orientation is defined in the literature as 
p m (z) = (pi) /[p 2 Pd(z)} : where (p 2 ) = / f±p d (z, n)p 2 z . 
The function p m was studied in Ref. [IT] for multipolar 
ions and it was found invariably below the free dipole 
value 1/3 in the SC limit (i.e. dipolar alignment parallel 
to the wall) and above this value in the WC limit (align- 
ment along the electrostatic field). We show in Fig. [3] 
that for a neutral interface, one has p m {z) < 1/3 as in 
the SC limit, i.e. the solvent molecules exhibit a ten- 
dency to align parallel to the wall over a distance f=a 2 A, 
that is, until image dipole forces vanish. We now define 
an effective dielectric permittivity function of the form 
e e ff(z) = l + 4:ir£BPoPd(z)/3 that solely accounts for the 
dipolar depletion. The comparison of i e ff(z) in Fig. [TJa 
with e(z) shows that the main effect of the dipolar align- 
ment close to the interface is a slight reduction of the 
local dielectric permittivity. However, it is seen that this 
effect is largely dominated by the solvent depletion. 

In the presence of a finite surface charge, Fig. pi) shows 
that interestingly, the function p m (z) exhibits a non- 



IV. CONCLUSIONS 
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(a) 

FIG. 3: (Color online) (a) Dipolar orientation profile eval- 
uated for various a s and the same model parameters as in 
Fig. n\ The black dashed reference line marks the freely ro- 
tating dipole case fJ, m (z) = 1/3. 



monotonous behavior. Namely, in the close vicinity of the 
surface, strong image dipole effects lead to a net dipolar 
alignment parallel to the wall. However, above a char- 
acteristic distance from the dielectric interface where the 
surface charge induced electric field pq4>q(z) dominates 
the image dipole potential, n{z) exceeds 1/3 and reaches 
a peak where the dipoles exhibit the maximum tendency 
to align in the direction of the field, i.e. perpendicular 
to the dielectric wall. With increasing distance, one no- 
tices a reversal of this behavior where the image dipole 
potential dominates for a second time the electrostatic 
field. As a result, the solvent molecules exhibit again 
some tendency to align again parallel to the interface, 
but this regime gradually disappears with increasing sur- 
face charge. 



We have considered the dielectric discontinuity ef- 
fects on the differential capacitance of low dielectric sub- 
strates. To this aim, we derived an extended DPB equa- 
tion that can explicitly account for the interactions of 
solvent molecules with their electrostatic image. Within 
this approach, we showed that the overestimation of the 
experimental data by the GC capacitance is due to the 
inability of the latter to account for the solvent depletion 
effect driven by image dipole interactions. The prediction 
of the EDPB equation for the differential capacitance of 
monovalent electrolytes was compared with experimental 
data and good agreement was found. 

The EDPB formalism is a first order theoretical ap- 
proach in the explicit modeling of solvent molecules be- 
yond the MF level, and it has its limitations. Excluded 
volume [IB] and non-local dielectric effects that lead in 
MD simulations to intcrfacial structure formation [9] are 
absent. The theory could be extended by using more gen- 
eral trial kernels, but the analytical solution of the Debye- 
Huckel equation with a local dielectric permittivity and 
screening parameter is still an open problem. Further- 
more, the present formalism does not include multipolar 
moments, which are known to enhance the interfacial di- 
electric exclusion 1 17]. Hence, multipolar contributions 
are expected to further lower the capacitance curves in 
Fig.0 
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